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Stability chart of the triangular points in the elliptic 
restricted problem of three bodies 
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ABSTRACT 

The possible observations of Trojan-like extrasolar planets stimulate the deeper under- 
standing of the stability behaviour of the co-orbital resonant motion. By using Hill's 
equations and the energy-rate method an analysis of the stability map of the elliptic 
restricted three-body problem is performed. Regions of the /i — e parameter plane are 
described numerically and related to the resonant frequencies of librational motion. 
Stability and instability can therefore be obtained by analysing the two independent 
frequency modes depending on system parameters. The key role of the long period 
libration in determining the structure of the stability is demonstrated and also a sta- 
bility mechanism is found that is responsible for extended life time of the test particle 
in the unstable domain of the stability map. 

Key words: celestial mechanics - planets and satellites: dynamical evolution and 
stability - methods: numerical, analytical 



I INTRODUCTION 

The stability of the Lagrangian triangular points 1/4 and L5 
in the elliptic restricted three-body problem (ERTBP) is an 
extensively studied topic of celestial mechanics. Recent di- 
rections in extrasolar planetary dynamics require to rethink 
our knowledge about the historical results of pure dynam- 
ical astronomy. The stability of motion around I/4 and L5 
depends on two parameters in the ERTBP, namely on the 
orbital eccentricity (e) and the mass parameter (/i) of the 
primaries. 

The ma in results on t he stabil i ty of Z/4 a nd I/5 were 
publis hed in lDanbvl l| 19641 ). iBemiettI (l965l ). and lTschauiie^ 
(|l97ll ). Nu merical studies about th e st ability was car- 
ried o ut bv lLohinger fc Dvo rak' (1993') and 'Markel los et al.l 

II 19951). It was sh own by Erdi ct al.. (2007 , 20091 ) and 
iRainai et all (|2010l ) that there are secondary resonances re- 
sponsi ble for the fin e stru cture of the stability chart. Re- 
cently [Sc^warz^^et^al] (|2012l ) demonstrated the stability phe- 
nomena in the spatial problem, too. 

In this work we present a fast and accurate method to 
construct the detailed stability map of the problem using 
Hill's equations. The results reveal the complex structure of 
the p — e plane and give evidence for peaks of longer escape 
times in the unstable region of the parameter space. 
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2 EQUATIONS OF MOTION AND METHODS 

The behaviour of the test particle near the L4 and L5 
points in the planar elliptic restricted three-body problem 
can bo described by the well-known equations of variation 
(Szebchclv, 19671) 



X — 2y — rcix, 
y" + 2x' — rC2y 



(la) 
(lb) 



where 



r = (1 -I- e cos v) 
m2 



mi + 7712 

g = 3/^(1 - n), 
= 1.5(1 + (-1)'^!^ 



(i = l,2) 



In Eqs. H]) the prime denotes the derivative with respect 
to the true anomaly v. The orbital parameter e and the 
dimensionless mass parameter jj, are involved in the problem, 
mi and m2 are the masses of the primaries. 



2.1 Hill's equations 

By a suitable transformation (|Tschauneilll97ll ) . the 4-th or- 
der system ^ can be separated into two independent 2-nd 
order systems. From these we can derive an equivalent set 
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of Hill's equations (|Meirelll980l . Il982l ): 

5i + Ji{v; e, = 0, 
(.2 + J2{v;e,y)£,2 



0, 



(2a) 
(2b) 



where Ji{v;^,e) are 27r-periodic functions containing the 
system parameters. Therefore the system ([2]) can be thought 
as of two decoupled harmonic oscillators with variable fre- 
quencies. The characteristic units in Eqs. ^ and ([2]) were 
chosen so that the primaries' period is 2n, i.e. the mean mo- 
tion is 1. For a detailed description of Eq. ^ see Appendix 
A. 



2.2 Energy-rate method 

The energy-rate method (ERM) (|jazaj|2003 ') gives the op- 
portunity, beside the well-known methods (Floquet theory, 
Poincare-Linstedt series) , to investigate the stability bound- 
ary of a system 



where 



5(0,0,t) = 



(3) 



(4) 



and f{x) is a single variable function, and g{x,x, t) is a non- 
linear periodic function of the time t. The functions / and g 
may depend on some parameters of the system. It is trivial 
from Appendix A that Eqs. ((2)1 can be written in the form 
of Eq. @ 



iz ~ - Ji{v\ ti, e)6 = 0, 



where 



Ji{v; fi, e) = rci 



3Qi -I- C2 

{') 

912 



+ 2. 



(i = l,2) 



(i) 
'722 



,'7i2 



(5) 



(6) 



Ji corresponds to excitation g(x, x, t) in Eq. ((Sj as one may 
check in Eqs. (lAlll . 

As is well-known if any non-conservative force is present 
in the system, the mechanical energy does not remain con- 
stant. In other words, due to the effect of g{x, x, t) in Eq. ([3]) 
for a given set of parameters, energy is either pumpin g into 
or w ithdrawing from the system. It can be verified Jjazail 
[2003) that in a system such as ((3| the energy integral has 
the following form 



E=f^iE)^fiix^+ lfix)xdt 



-xg(x,x,t) (7) 



where T{x) = i^/2 is the kinetic and V(x) — J f{x)dx is 
the potential energy for free motion. 

It is also known that the time derivative of the mechani- 
cal energy must be zero over one period for conservative and 
autonomous systems in a steady state periodic cycle. In or- 
der to find the parameters that belong to stable or unstable 
solutions we may evaluate the energy over a period 



-Bav — rp 



Edt 



(8) 



where T is the period of the external force g. If Bav > 0, 
provided that / is a restoring force, then energy is inserted 
into the system, and Eq. ([3}) is unstable. But, if _Eav < 
0, energy is subtracted and the system is stable. On the 
common boundary _Eav = 0. 
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Figure 1. Stability map in the fj, — e plane. Top: different fea- 
tures of motion around L4 and L5 governed by Eq. I|2ap . Bottom: 
stability domains corresponding to Eq. l l2bll . The contour bar in- 
dicates the magnitude of the average energy for given pairs of the 
parameters (/J, e). The integration time is 2it, one period of the 
primaries. The number of the parameter pairs is set to 10^, i.e. 
1000 in each direction with equal step size. (Colour figures are 
available online.) 



3 STABILITY ANALYSIS 

In the following the stability chart of the L4 and L5 points 
is studied using the ERM. A common way in stability anal- 
ysis is to identify the characteristic roots of the problem 
that characterize the dynamically different domains of the 
parameter plane. The ERM allows us, if the equations are in 
a suitable format, a simple and fast method to distinguish 
the sets of parameters corresponding to stable or unstable 
dynamics. 

In order to get the classical stability map of the /i — 
e plane one has to compute Eq. (|8} numerically for both 
Eqs. (O 



(i = l,2). 



(9) 



Note that equations ([2]) have a repulsive force, /(^) — 
—2^, (see ((H) responsible for free motion. Therefore the sta- 
ble solution for some (/i, e) pairs can be realized when the 
averaged energy _Eav > 0. Consequently, unstable motion 
appears in the opposite case, iJav < 0. 

Since Eqs. ((2)1 are linear ODEs with periodic coeffi- 
cients, the choice of initial conditions is arbitrary, i.e. they 
do not affect the stability of the solution. Throughout this 
work initial conditions (Ci,?0 = (1,0) are used. The control 
parameter that measures the accuracy iSav = 0, is set to 
e = 5 - 10"^ 

Figure [Ha) depicts the classical stability map of the 
ERTBP obtained from Eq. ^ by the ERM. The light gray 
region (orange online), where i5av > for one period, cor- 
responds to librational motion around L4 and L5. The dark 
gray area (blue) represents the negative energy domain re- 
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lated to the set of parameters of unstable motion. The black 
part of the figure (dark blue), situated at the bottom right 
quarter, belongs to parameters where Eqs. ((2| have no so- 
lutions. For better visualization (log-scale) we cut the pos- 
itive energies at 10^ and set the negative values to 10"^", 
additionally, where the problem has no solution to 10~^*. 
Therefore, the contour bar is organized as follows: the light 
gray colour (orange, 10"^" -Eav ^ 10^) represents libra- 
tion, the dark gray (blue, iSav = 10~^°) denotes unstable 
motion, and the black domain (dark blue, -Eav = 10~^*) is 
for complex frequencies. The border of the stabilit y domain 
is at _Eav ~ 0. The charecteristic roots of Eqs. dBennettl 
Il965l) devide the stability chart exactly into the same parts. 

Figure [TJb) shows the stability chart provided by 
Eq. (j2b[) . It can be seen that the librational domain (light 
gray (orange online)) is more extended than in panel (a), 
however, the bottom right part (black/dark blue) is exactly 
the same. If we want to identify the stability region around 
the Li and I/5 points in the ERTBP determined by Hill's 
equations, we need the intersection of all the individual sta- 
ble regions in Fig. [TJa) and[TJb). One can immediately verify 
that the upper panel encompasses the domain in question. 
In other words, this shows that merely the first Hill's equa- 
tion in Eqs. ([2| is responsible for the well-known structure 
of the /i — e plane. The possible reason will be discussed in 
Section 5. Accordingly, we can state that the stability chart 
obtained from Eq. ([2]) is consistent with earlier analysis. 

One of many advantages of the ERM is that if we in- 
tegrate the equation(s) of motion until t — nT where n 
is an integer number and T is the period of the external 
force, new, so-called "splitting" curves appear on the stabil- 
ity map. (The name comes from the fact that these resonant 
periodic lines split the stable regions into several parts.) Ac- 
tually, these curves are super- and sub-harmonic oscillations 
and describe higher order periodic responses of the paramet- 
ric equations. 

Figure [5] shows the average energy levels at t = Stt, 
i.e. for 4 revolutions of the primaries. The similarity with 
Fig. Ul is obvious but we have several new splitting curves 
inside the stable regions, both in panel (a) and (b). In Sec- 
tion 4 it will be shown that these resonances cor respond to 
the ty pe A and D secondary resonances defined in lErdi et al.l 
(|2009l ). 



4 RESONANT FREQUENCIES 

It was shown bv lErdi et all (|2009l 'l that the fine structure of 
the /Lt— e plane is organized by secondary resonances between 
the four frequencies in the ERTBP. Two of them, the short 
Us and long n; frequencies (corresponding to the short and 
long periodic librations around L4 and L5), are present also 
in the circular problem, the additional two (in normalized 
units: 1 — Jis, 1 — nj) appear when the primaries revolve on 
elliptic orbits. 

In general case. Hill's equations ((2]) can be considered as 
of two independent harmonic oscillators with periodic coeffi- 
cients, i.e. with variable frequencies. Eqs ([2]) describe the two 
fundamental frequencies in ERTBP, namely the frequencies 
corresponding to the short and long periodic components of 
the libration around L4 and L5. When the problem is circu- 
lar, Eqs. ([2]) become simple harmonic oscillators and serve 
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Figure 2. The /i — e plane for integration time t = Sir. Top and 
bottom: stability zones derived from Eqs. I|2a|l and l)2b| l. respec- 
tively. Increasing the computation time by the multiplication of 
the excitation period, resonant curves turn up inside the stable 
regions. 



the well-known frequency diagram of the problem. (See Ap- 
pendix A, Fig. lAin 

The ERM apphed to the ERTBP can be used to identify 
the secondary resonances between the librational frequencies 
of the test particle and the primaries' motion. These are 
the resonant curves on the stability c hart de fined as A-type 
(l — ni) : Hi and D-type Ug : (1 — n^) in lErdi et al.. (|200Sjh . In- 
vestigating the problem through Eqs. ^ one do not need to 
involve the two "additional" normalized frequencies coming 
from elliptic ity. Let us consider a quick example. The All 
resonance inl Erdi et al.l (j2009') corresponds to the nearly V- 
shape boundary line of the stability zone. The notation All 
comes from (1 — n;):n; = l:l. However, after a simple rear- 
rangement we get 1/2 = ni, i.e. the frequency is equal to the 
half of the normalized frequency of the primaries. In other 
words, the long period of the librating infinitesimal parti- 
cle alon g the All res onance curve is twice as the primaries' 
period l|Danbvlll963 ). Moreover, this rule holds for any A 
and D-type resonance. If we integrate Eqs. ^ for t = 2-n 
and plot the sets of parameters corresponding to i5av = 
condition, the stability map with secondary resonances, of 
types A and D, can be drawn from. 

In Figure [3] the ji — e plane is depicted for various in- 
tegration times. Panel (a) contains the "original" stability 
chart of librational motion in the ERTBP with agreement 
of the energy plot of Fig. [T] Dark blue curves (color online, 
hereafter dark curves) correspond to the long period motion 
(ni) while the light blue (color online, hereafter light curves) 
represent the short period (n^). We note again, these res- 
onant lines indicate the two components of motion around 
the points L4 and L5 . Accordingly, the nearly V-shaped dark 
curves describe the 2:1 commensurability of the test parti- 
cle's motion, in direction ^1, with the primaries' period. The 
light line represents the resonant component in the other di- 
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Table 1. A- type resonances in Figure (Sfd). 



Resonance starting point panel (s) 

8:1 0.00228 

4:1 0.00876 b 

8:3 0.01823 

2:1 0.02859 a,b,c 

8:5 0.03660 



Second column contains the mass parameters along the e = axis 
where the resonant lines originate. Third column gives the other 
panels where the same resonance appears. 

rection ^2- Along these curves the average energy _Eav, during 
one period is zero. This picture itself does not say anything 
about the stability of the two-dimensional motion but only 
the resonant frequencies. We also need the characteristics of 
Figs. [1] and [2] to decide whether the motion is stable or not. 
Fig. [3Ja) tells us only the border of the stability region. It 
has already been shown in Section 3 that libration can occur 
under the dark curve, where both components of Eqs. ((2]) 
are stable. 

Other secondary resonances show up on the stability 
map when the integration time is doubled, as in Fig. El^b). 
In this panel two extra lines appear, beside the well-known 
stability border, corresponding to the 4:1 mean motion res- 
onances . The plot was obtained as follows: the integra- 
tion time was 47r, two periods of the primaries, the (/i, e) 
pairs were stored when the energy rate was zero during 
this time. Now, the harmonic resonant curves are the new 
lines and the old borders are the sub-harmomcs of them 
(2:1 res.). This behaviour is defined as splitting in the 
ERM, where resonant curves split the stability regions into 
two parts ((jaza r 2004). The 4:1 resonant curve originating 
from (fi, e) = (0.0 087, 0.0) corresponds to A31 resonance in 
lErdi et all (|2009l ). To demonstrate the feature of the reso- 
nances Fig. |4] shows two individual solutions of A- type res- 
onances, 2:1 and 4:1. 

If one continues the procedure, more and more splitting 
curves turn up on the stability map. In panel (c) {tint ~ 67r, 
3 revolutions of the primaries) the fund amental res o nance s 
are 6:1 - {^^,e) = (0.00397,0.0) (A51 in lErdi et al] (|2009l )) 
and (fj., e) — (0.0369, 0.0) for type D resonance. The sub- 
harmonics describe the 3:1 (A21) and 2:1 (All) resonances. 

In addition two interesting curves appear in panel (c). 
Curve a originates from {fi, e) = (0.0381, 0.0), curve d comes 
from the stability region of ^2 and terminates where curve 
a ends, about at (p, e) = (0.042, 0.202). Integrating Eqs. ^ 
along these curves, one can find that these librations corre- 
spond to the super-harmonics of the resonances (6:1). These, 
in question, are the 2:3 commensurabilities. The role of the 
meeting point is discussed in Appendix B. 

Figure [3jd) demonstrates the resonant frequencies for 
tint = 4 revolutions of the primaries (Stt). As one may ex- 
pect various new resonances appeared on the ^ — e plane. 
Table 1 collects the resonances corresponding to dark curves 
on panel (d). D-type resonances are light curves of similar 
types. 

To show how complex the stability chart can be, the 
fj. — e plane for tint ~ 207r is plotted in Figure [5] The 10:1 
resonances and all their sub- and super-harmonics are de- 
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mass parameter 

Figure 3. Resonant frequencies at different integration times. 
Computations were performed from top to bottom for times: 
27r,47r, Ott, and 8tt. In panel (b)-(d) splitting curves appear. For 
detailed explanation see the text. 

picted. An even more complete and dense structure could 
be obtained by overplotting the charts corresponding to dif- 
ferent integration times. 



5 SUMMARY AND CONCLUSIONS 

In the present work the stability analysis of the ERTBP was 
performed. Using Hill's equations we showed that the prob- 
lem can be described as two independent harmonic oscilla- 
tors with time dependent frequencies. It was also demon- 
strated that using this approach, combined with the ERM, 
one can find the same stability chart as many authors pub- 
lished previously. 

Beside many advantages, the energy-rate method allows 
us to obtain the resonant curves on the fi—e plane with arbi- 
trary accuracy e (obviously depending on numerical trunca- 
tion errors). Moreover, it is not restricted to moderate values 
of the parameters as common perturbation technics. 

Investigating the librational motion around the L4 and 
L5 points in the ERTBP by two independent Hill equations, 
two unresolved problems in stability analysis can be inter- 
preted 

(i) Figs. [T] and [5] indicate that the long period component 
(frequencies from Eq. (|2a|l ) of libration is responsible for the 
global feature of stability. Indeed, the A-type resonances are 



Stability of the triangular Lagrangian points 5 




10 20 30 40 50 
time 



Figure 4. Two individual solutions of Eq. Il2ap . The dashed 
line corresponds to the 2:1 resonant frequency for the parame- 
ter pair {fi,e) = (0.02312,0.1), see the circle on the dark curve 
in Fig.[3lb). While the slolid curve denotes the 4:1 resonance for 
(/i, e) = (0.00838, 0.1), this latter point is marked with a cross in 
the stability chart, see also Fig.[3lb). The dashed and dotted line 
represents t=25.12 (Sir), i.e. 4 revolutions of the primaries. 
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Figure 5. Stability chart for t = 207r. This plot contains the 10:1 
resonant curves and all sub- and super-harmonics of them. 

situated in the intersection part of the two stable regions 
(long and short period librations), they are therefore the 
backb one of the two-dime nsional stable motion. 

(ii) iRainai et all (|2010l ) showed that there are significant 
peaks in the unstable domain responsible for longer escape 
times. In other words, there may exist some processes that 
stretches the lifetime of t he test partic l e arou nd L4 and L5. 
Comparing their result (|Rainai et al.l l|2010l ). Fig. 1) with 
Figure [S] one can see that the longer escape times appear 
exactly where the resonant frequency curves, corresponding 
to the shorter frequencies, are situated. An explanation of 
this behaviour, based on the present work, can be that the 
resonant curves in the stable domain of ^2, that are actually 
invisible on the classical fi — e stability chart, may partially 
stabilize the whole librational motion around the triangular 
Lagrangian points. 

A direct consequence of this study can be, in future, to 
extend the investigation to the spatial problem, considering 
the third equation of motion. Interestingly, the third vari- 
ational equation beside the x and y components of ^ is 
already in form of Hill's equation, therefore, only a simple 
variable change is needed to study the extended ERTBP. 

The reader may ask why are the B-type resonances 



(|Erdi et al.ll2009l ) omitted from the present study. Since the 
ERM cannot be used to identify the exact solution of Hill's 
equations and therefore one cannot compare the indi- 
vidual frequencies, Us : ni (perhaps the most important) 
resonances were not examined in this work thoroughly. How- 
ever, an interesting analytical approach of Bll resonance is 
presented in Appendix B. 
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APPENDIX A: HILL EQUATIONS IN THE 
RTBP 

Following Meire 's (lMeirelll982D derivation, the final form of 
Hill's equations in the restricted three-body problem are 

+ Mv;e,^i)^, = 0, i={l,2) (Al) 
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where 



with 



Ci ■- 
9 ■- 



3Qi + C2 , „ 
'■ci + 2 (-)— + 3 

<ll2 



922 



912 



1.5 [i + {-iy^T^ 



3/i(l - fi), 

1712 
mi + 7712 ' 



Qi = 0.5 [l + (-l)V + 3ecosuj 



k 

(1) (2) '"C , 1, 2 o 

=912 — ^ = '^2 + e COS V — -ke cos 2«, 



§22^ ~ — iesinu(l — fcecosu) 



0) 

u 

0) 
4) 
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Figure Al. Librational frequencies in the circular RTBP. The 
darli blue (online) curve corresponds to the long period solution 
uii or n;. The light blue (online) curve depicts the short period 
component UJ2 or ng. Eqs. | |A2| | provide the exact frequencies 
depending only on the mass parameter as e = 0. 



02 = -(2c2 + l-c). 

In Eqs. jAl^ e denotes the eccentricity, v is the true anomaly, 
mi and m2 are the masses of the primaries, and r — (1 + 
ecosv)~^ is the mutual distance between them. 

One can obtain the motion around L4 and L5 as a com- 
position of two harmonic oscillations governed by Eqs. (|A1[) . 
Obviously, the form of equations shows that the frequencies 
of the motions are y/Ji. However, the solution is trivial only 
in the case when the frequencies do not depend on time, i.e. 
the orbital parameter is zero. After some computations, one 
can easily verify that in this case (e = 0) the frequencies 



(A2) 



where 



AM 



2.^[f 



3 + 1(^/1^+ (-1)"^/!^) 



i + |^r^+(-i)4^r^- 

Figure IAT] demonstrates oji and u}2 vs. the mass param- 
eter /i. It can be seen that the curves exactly match the 
well-known classical results of the circular RTBP. The lower 
frequency uji represents the long period solution while UJ2 
correspondes to short time oscillations. The latter tends to 
the frequency of the primaries (here normalized to 1) in the 
limit case when ^ — ^ 0. 



APPENDIX B: THE Bll RESONANCE 

In the elliptic case the frequencies of libration Eqs. HAlf) are 
not so obvious. In this case we have second order linear dif- 
ferential equations with time dependent coefficients. For an 
arbitrary given periodic coefficient, there are no general pro- 
cedures to express the characteristic exponents, unless the 
funda mental matrbc of the proble m is obtained by quadra- 
tures (|Arnoldll 19781 : lYoshid al ll984l ). Therefore, it seems to be 
a hard task to find any analytical expression for the period, 
and therefore the frequency Ui, of the solution ^i{v). 
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Figure Bl. The rij, : n; = 1 : 1 frequency curve (Bll) as the 
border of stability. A simple averaging by using Eqs. IIBll l gives 
the exact curve of the Bll resonance. This curve indicates the 
parameter pairs {fi, e) where lJ\ = tUJ- The resonant curves All, 
Dll, and Bll meet at one point P. Additionally, several resonant 
points can be identified for different mass parameters at e = 0. 
Crosses along the /i-axis denote the following B type resonances, 
from left to right: 5:1, 4:1, 3:1, 2:1, 3:2, 4:3, 1:1, respectively. 



Let us consider, however, a naive analytical approxima- 
tion to find the resonant curves corresponding to the B-type 
resonances. The B-type resonances are defined as: ria : nj, 
where and n; denote the short and long period frequen- 
cies, respectively lErdi et al.ll2009h . These frequencies, that 
are actually varying with time, can be easily obtained by 
evaluating the Ji coefficients in Eqs. (|A1[) . Suppose one can 
describe the solution with frequency uji changing around the 
average frequency ZUI 



dv : 



-J2VMh^^n^ (Bl) 



TV 



where T is the period of Ji, /i is a stepsize, and A*' is the 
number of terms during one period. Consider these averaged 
frequencies the skeleton of individual librations around L4 
and L5. In this sense, B-type resonances can be evaluated 
simply analytically by using the formulae given by Eqs. (|B1|) . 
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Figure B2. Various A and D-type resonances near to Bll reso- 
nance. The panels (a)-(c) refer to different integration times, 67r, 
Stt, IOtt. (Magnified parts of Figs.|3lc),|3ld), and[5l respectively.) 
On each panel the touching point is marked with a circle. At these 
points the two librational frequencies n,, and n; are equal. Panel 
(d) verifies that the Bll resonance obtained by using Eqs. lIBlll 
contains the circles. 



Just for the interesting coincidence with Erdi's work, 
who used Fourier analysis to find the exact period, the Bll, 
i.e. Us : ni = \ : \, resonance is plotted in Figure [Bl] As is 
well-seen the curve originates at (/i, e) = (0.0385, 0.0) and 
follows perfectly the border of stability, touching the critical 
point (/i, e) = (0.0414, 0.314) and draw out the border line, 
seen also in Fig. [T] that separates two different kinds of 
solutions (Eqs. (lAlfl ). No analytical solution in li near theory 
was f ound so far that can describe this behaviour jSzebehelvl 
[1963). 

At this point we should explain the role of the pair of 
the All and Dll (P in Fig. IB1|) resonant curves that touch 
each other in a common point. These A and D-type reso- 
nances describe a well-defined secondary resonance between 
the librating particle and the primaries' motion. However, at 
that particular point P{fi,e), where they meet, the two in- 
dividual libration frequencies become equal. Therefore, the 
motion at this point also satisfies the condition of Bll res- 
onance, i.e. uji = UJ2. For different integration times new 
A and D-type resonant curves appear on the fi — e plane 
that meet in one common point, see Fig. IB2I These com- 
mon points are necessarily situated along the Bll resonance 
curve due to the condition uji = cj2. 

Finally, we also note that this method provides the right 
places of the secondary resonances for the circular RTBP. 
That i s, all the resona nt frequencies (with A,B,C,D,E,F- 
types, (|Erdi et all boogP l can be identified by the analyt- 
ical computation Eqs. (|A2[) along the /^-axis, see Fig. IBll 
The reason of this is that the coefficients in Eqs. (|A1|) no 
longer depend on time in the case e = 0, consequently Ji 
describe the exactly constant frequencies of the librational 
components, since in the circular case Ji ~ ijjf. 
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